Collective motion and nonequilibrium cluster formation in colonies of gliding bacteria 
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We characterize cell motion in experiments and show that the transition to collective motion 
in colonies of gliding bacterial cells confined to a monolayer appears through the organization of 
cells into larger moving clusters. Collective motion by nonequilibrium cluster formation is detected 
for a critical cell packing fraction around 17%. This transition is characterized by a scale-free 
power-law cluster size distribution, with an exponent 0.88 ± 0.07, and the appearance of giant 
number fluctuations. Our flndings are in quantitative agreement with simulations of self-propelled 
rods. This suggests that the interplay of self-propulsion of bacteria and the rod-shape of bacteria is 
sufBcient to induce collective motion. 
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In many microorganisms, the transition from single cell 
to multicellular behavior involves the onset of collective 
motion, that is characterized by the formation of large 
cell clusters that move in a coordinated manner. An 
open question is by which mechanisms such cellular orga- 
nization is achieved. There are several examples of coor- 
dinated cell motion resulting from intercellular signaling 
systems. In the developmental cycle of Dictyostelium dis- 
coideum, a diffusive chemoattractant guides cell motion 
and leads to complex pattern formation [l|, [2| , while in 
Myxococcus xanthus a signaling system, that requires cell- 
to-cell contact, coordinates cell movements and gives rise 
to rippling patterns Q • In absence of a signaling system, 
spatial cellular organization can result from density de- 
pendent difFusivity, as suggested to occur in Escherichia 
coll and Salmonella typhimirium In other cases, 

large-scale coherent patterns are believed to emerge from 
hydrodynamic interactions, as in swimming bacteria like 
Bacillus suhtilis [3]. Thus, collective motion of cellular 
populations typically involves physical and biochemical 
interactions between cells |6||. A purely physical mech- 
anism has been recently analyzed in simulations of self- 
propelled rods 0, Q J ill which the interplay of active rod 
motion and steric interactions due to volume exclusion 
leads to the formation of moving clusters for sufficiently 
high densities. Since such collective movement was found 
to be absent in the equivalent equilibrium system of diffu- 
sive rods, this phenomenon can be considered as nonequi- 
librium cluster formation. Experiments with granular 
particles, i.e. artificial self-propelled rods, con&med that 
such a physical mechanism is indeed enough to produce a 
variety of collective motion pattern Q . Since many bac- 
teria are self-propelled and have rod cell shape, simple 



steric interactions may be sufficient to induce collective 
motion even without the additional impact of biochemi- 
cal signals. The aim of this study was to test this idea in 
experiments using M. xanthus as model system. 

M. xanthus is a gliding bacterium that has been re- 
peatedlyused to study pattern formation [l3| , social be- 
havior ll| , and motility 12[ . Locomotion of M. xanthus 
involves two different motility systems: the S-motility 
system, which depends on type IV pili and requires cell- 
to-cell contact [13| and the A-motility system that al- 
lows individual cells to move [l^. Force generation by 
the A-system has been suggested to rely on either slime 
secretion from the lagging pole [l3| , or on focal adhesion 
complexes [3]. Cells occasionally reverse their gliding 
direction and the reversal frequency is controlled by the 
frz chemosensory system [l^ [3| • 

Here, we report on the bacterial self-organization in 
experiments with a M. xanthus mutant (SA2407) which 
only moves by means of the A-motility system and is vir- 
tually unable to reverse. Complex interactions like social 
motility mediated by pili and cellular reversal are absent 
in this mutant. Hence, the experiments are suitable to 
test the theoretical hypothesis that collective motion can 
emerge from the combination of active motion of the cells 
and steric interactions due to volume exclusion 0, S] . 

Experiments consisted in spotting a drop of the desired 
cell density on an agar surface, and subsequently the bac- 
terial colony evolution was monitored by taking images 
of cell arrangements every 30 min for a total of 8 hrs. 
To follow the detailed dynamics of cell arrangements, we 
also made time-lapse recordings of about 30 min, with 
successive frames taken every 30 seconds. Control exper- 
iments showed that SA2407 cells do not reverse for the 
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FIG. 1: Clustering of SA2407 cells, (a): Statistics of cell 
length-to-width aspect ratio and dispersion of the data in the 
length- width plane (inset), (b): Cells align their orientations 
upon collision, (c), (d): Myxobacterial cells form moving 
clusters. Arrows indicate the direction of motion of the mov- 
ing clusters; the time interval between snapshots is 15 min 
{rj = 0.11). (e)-(h): The dynamical clustering process reaches 
a steady state that strongly depends on cell density. Typical 
snapshots corresponding to packing fractions rj — 0.06 in (e) 
0.1 in (f), 0.16 in (g) and 0.24 in (h). 



time scale of our experiments, whereas the isogenic frz+ 
strain reversed with a mean reversal period of ~ 10 min. 
We found that individual cells glide at an average speed of 
V = 3.10 ±0.35 ^m/min, and exhibit an average width of 
about W = 0.7 /J,m and an average length of L = 6.3 fim, 
resulting in a mean aspect ratio of k = L/W = 8.9 ± 1.95 
(Fig. [Ija)) and a cell covering an average area a — 4.4 
firri^. Experiments were confined to packing fractions 
smaller than 0.26, where the cell dynamics is restricted 
to a monolayer - at larger packing fractions percolation of 
clusters as well as formation of multilayers are observed. 
The packing fraction i] is defined as -q = pa, with p the 
(two-dimensional) cell density and a the average covering 



area of a bacterium given above. 

We found that under these conditions over time cells 
organized into moving clusters. Time-lapse record- 
ings showed that collision of cells leads to alignment 
(Fig. mb)). When the interaction is such that cells end 
up parallel to each other and move in the same direction, 
they migrate together for a long time (typically > 15 
min) 19]. Eventually, successive collisions allow a small 
initial cluster to grow in size. In the individual clusters, 
cells are aligned in parallel to each other and arranged in 
a head-to-tail manner, as previously described [2^. In a 
cluster, cells move in the same direction. Cluster-cluster 
collision typically leads to cluster fusion (Fig. [IIc),(d)), 
whereas splitting and break-up of clusters rarely occur. 
On the other hand, individual cells on the border of 
a cluster often spontaneously escape from the cluster. 
These two effects, cluster growth due to cluster-cluster 
collision and cluster shrinkage, mainly due to cells escap- 
ing from the cluster boundary, compete and give rise to a 
non-equilibrium cluster size distribution (CSD). Typical 
snapshots of cell arrangements for various packing frac- 
tion rj at the steady state are shown in Figs.[lje)-(h) and 
reveal a strong increase of the cluster size for increasing 
packing fraction rj. 

The CSD - p(m, t) - indicates the probability of a bac- 
terium to be in a cluster of size m at time t. Note that 
along the text, the term CSD always refers to this defini- 
tion. The cluster size distribution can be alternatively 
defined as the number rimit) of clusters of size m at 
time t. There is a simple relation between these two 
definitions: Pm{t) oc mnm{t). In experiments we have 
observed that the CSD mainly depends on the packing 
fraction 77. Hence, for all snapshots first the packing frac- 
tion was determined. Then, images with similar packing 
fraction rj were compared and the CSD was reconstructed 
by determing the CSD for all images within a finite in- 
terval of the packing fraction. The CSD p{m,t) reaches 
a steady state p{m) after some transient time. The du- 
ration of this transient depends on the packing fraction 
rj and is below 120 min for all rj < 0.2, see [l^. This in- 
dicates that the (non-equilibrium in the thermodynami- 
cal sense) clustering process evolves towards a dynamical 
equilibrium, where the process of formation of cell clus- 
ters of a given size is balanced by events in which clusters 
of this size disappear by either fusing with other clusters 
or by loosing individual cells from their boundary. This 
cluster dynamics is in sharp contrast with cell cluster 
formation driven by differential cell adhesion and/or cell 
proliferation as observed, for instance, in cancer cell ex- 
periments 21], where the CSD never reaches a steady 
state. 

The steady-state CSD p{m) strongly depends on the 
packing fraction 77, with more and more cells moving in 
larger clusters for increasing packing fraction rj. This is 
evident in Fig. [51 where we observe that at small values 
of 77, p[m) exhibits a monotonic decay with m, while at 
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FIG. 2: Cluster statistics of SA2407 cells. The figure shows 
the steady-state cluster size distribution (CSD) p{m) for var- 
ious packing fractions 77. Notice the qualitative functional 
change exhibited by the CSD p(m). At the critical point 
rjc ~ 0.17, p{m) ~ m'"' with 70 = 0.88. We define this tran- 
sition point as the onset of collective motion, see text. 



large tj values, p(m) is non-monotonic, with an additional 
peak at large cluster sizes. The CSD here was deter- 
mined at a fixed time (450 minutes) after the beginning 
of each experiment; control measurements at other times 
(360 minutes, 540 minutes) revealed practically identical 
behavior. We interpret the presence of a peak at large 
values of m, observed for large ry values, as the emergence 
of collective motion resulting in formation of large clus- 
ters of bacteria moving in a coordinated fashion. This 
phase is characterized by the existence of large clusters 
that are reflected by the emergence of a local maximum 
in the CSD, see Fig. [21 The transition is evident by the 
functional change displayed by p{m), monotonically de- 
creasing with m for small values of 77, while exhibiting 
a local maximum at large rj values. At a critical value 
rjc = 0.17 ± 0.02 that separates different regimes of be- 
havior, the CSD can be approximated by p(m) cx m~'^°, 
with 7o = 0.88±0.07. In summary, for rj < rye, the scaling 
of p{m) takes the form: 



p(m) oc m exp(— m/mg) , 



(1) 



while for i] > rjc, the scaling is: 



p{m) (X m ^'^ exp(— m/r7ii) + Cm''^ exp(— 771/7712) , (2) 

with 71, 72, 7771, and C constants that depend on 77. 
In Eq. IT]), mo is a function of 77 and increases as 77c is 
approached from below. Thus, at the critical packing 
fraction 77c, p{m) can be approximated by a power-law 
as long as m is much smaller than the total number of 
cells N in the system. Eqs. ^ and ([2]) were obtained in 
self-propelled rod simulations ^, |3fl] and used here to fit 



the data in Fig. [2l Control experiments with non-motile 
cells do not exhibit power-law behavior in the CSD for 
larger packing fraction 22]. Hence, we conclude that 



without active motion of cells no comparable transition 
to collective motion occurs. In other words, active motion 
is required for the dynamical self-assembly of cells. 

It is interesting to observe that simulations with self- 
propelled rods PjliO] exhibit a very similar behavior of 
p{m). Moreover, the exponent 7n takes on similar values 
as in the experiment. Data in [3,|33| give 70 = 0.95±0.05, 
and Yang et al. Q report on similar simulations, obtain- 
ing values for 70 in the range from 0.95 to 1.35. On 
the other hand, the mean-field theory for the cluster-size 
distribution introduced in [Tj gives an exponent of 1.3, 
which is much larger than the experimental value mea- 
sured here. This theory can be extended to account for 
elongated rather than circular cluster shapes. In the ex- 
tended theory, predictions for the exponent 70 depend on 
the functional form of the coagulation and fragmentation 
kernel. Upon a series of assumption, the theory predicts, 
for elongated clusters, an exponent of 0.85 [31 1. 

We have also characterized the number fluctuations of 
SA2407 cells: (A77(/)) = (77(0^) - (n(0)^ where 77(/) 
denotes the number of cells in a box of linear size I. It 
can be shown that in general: 



(A77(0) cx (77(0)^ 



(3) 



with {n{l)) = pp. Thus, the quantity An is a mea- 
sure of the distribution of cells in space. Normal fluc- 
tuations correspond to /3 ~ 1/2, while for giant fluc- 
tuations /3 > 1/2. It has been argued that systems of 
self-propelled particles exhibit in their (orientational) or- 
dered phase giant number fluctuations, which are often 
considered a signature of non-equilibrium [23| . It has 
been shown recently that self-propelled particles with ap - 
olar alignment effectively exhibit such fluctuations [24 1. 
Fig. [3] shows that (A77) is a function of the density p, re- 
spectively, 77. At low values of 77, number fluctuations are 
consistent with normal fluctuations. We observe, never- 
theless, that for small 77, (Ati) exhibits a crossover from 
a regime characterized by an exponent close to 0.8 for 
small (77) to an asymptotic regime for large (n) charac- 
terized by an exponent 0.5 as expected for normal num- 
ber fluctuations. As 77 is increased towards 77c, cells ex- 
hibit giant number fluctuations though there is a lack of 
global orientational order. Morevover, for rj > rjc num- 
ber fluctuations are characterized by the same exponent 
/? = 0.8 ±0.05 (Fig. [3]). Interestingly, this value coincides 
with the exponent reported in [24|. These flndings sug- 
gests that in the experiments giant number fluctuations 
are connected to the transition observed in the cluster 
size statistics. 

We have reported on the bacterial cell self-organization 
in experiments with a M. xanthus mutant which only 
moves by means of the A-motility system and is vitually 
unable to reverse. We found that these bacteria exhibit a 
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FIG. 3: Number fluctuations of SA2407 cells - at low packing 
fractions rj, number fluctuations An are consistent to what 
is expected for normal fluctuations. For large values of ry, 
i.e., for rj > rjc, number fluctuations are giant, with a critical 
exponent 0.85 ± 0.05. 



transition at a critical packing fraction of rj = 0.17 char- 
acterized by the emergence of a power-law cluster size 
distribution, with an exponent 70 = 0.88 ± 0.07, and gi- 
ant number fluctuations, with exponent (3 — 0.8±0.05, in 
the absence of global orientational order. The observed 
change in the spatial organization of cells with increasing 
packing fractions resembles that obtained in simulations 
with excluded volume interactions among self-propelled 
rods 0. Moreover, we observe that there is a surprising 
similarity between the obtained statistics for the glid- 
ing M. xanthus mutant studied here and the recently 
reported results for the swimming bacterium B. suh- 
tilis [25'! . These observations raise the questions whether 
(i) the spatial organization at the onset of collective mo- 
tion may have universal features and (ii) if such features 
may indeed be linked to the simple physical paradigm 
of self-propelled rods. The latter speculation finds addi- 
tion tentative support by very recent experimental mea- 
surements on the nature and range of cell-cell interac- 
tions in suspensions of E. coli, where it was found that 
the emergence of large-scale patterns in bacterial sus- 
pensions is likely to be dominated by simple physical 
collisions among bacteria rather than by hydrodynami- 
cally induced long-range dipole-forces [26|. Ultimately, 
only further experimental studies of bacterial colonies 
and biofilms can reveal whether universal features are 
actually present at the onset of collective motion in a 
wide-range of cellular systems. 

In summary, we have shown that the cluster size distri- 
bution exhibits a qualitative (non-equilibrium) transition 
(reminiscent to a gelation transition [27'|) from an expo- 
nential decaying shape at small 77 to a power-law shape 
with an additional peak at large rj. We have suggested 
to use this characteristic transition point in the CSD as 
definition for the onset of collective motion. According 



to this definition, collective motion implies the forma- 
tion of large moving clusters, with cells sharing the same 
moving direction inside the clusters. Since in the exper- 
iments clusters do not exhibit a moving directional pref- 
erence, myxobacteria exhibit collective motion without 
global orientational ordering. The reported transition to 
collective motion (via clustering) is hence qualitatively 
different from the transition to global orientational order 
reported in the Vicsek model and its variants [IS, HE 113 ■ 
It would be interesting to explore the relation between 
both transitions, i.e., via clustering and global orien- 
tational order, in such minimal models where recently 
nonequilibrium cluster formation similar to the one ob- 
served here for colonies of gliding bacteria or for simula- 
tions of hard rods has been reported [s^, 32 1. Moreover, 



we have found that the number fluctuations in the ex- 
periments are normal in the limit of large numbers {n) 
for densities below the onset of collective motion. In con- 
trast, above the transition giant number fluctuations are 
found. The cluster-size statistics as well as the number 
fluctuation show both a distinct qualitative change at the 
onset of collective motion. Consequently, both measures 
are suitable to describe the onset of collective motion in 
large groups of microorganisms. 

Finally, the agreement between the cluster statistics 
obtained in the experiments and in earlier simulations of 
self-propelled hard rods 0, [l] , suggests that the interplay 
of active motion and volume exclusion is sufficient to ex- 
plain the collective behavior of the bacteria considered 
here. 
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